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We study numerically a parametrically driven discrete nonlinear Schrodinger equation modelling periodically 
curved waveguide arrays. We show that discrete surface solitons persist, but their threshold power is altered by 
the drive. There are critical drives at which the threshold values vanish. We also show that parametric drives 
can create resonance with a phonon making a new barrier for discrete solitons. By calculating the corresponding 
Floquet multipliers, wc find that the stability of symmetric and antisymmetric off-side discrete surface solitons 
switches approximately at the critical drives for thresholdless solitons. © 2013 Optical Society of America 
OCIS codes: 030.1640, 190.4420 



Parametric drives have been proposed as a means to 
control localized waves in a linear Schrodinger system, 
which in the undriven case would simply disperse [1]. 
Such drives are created among others in periodically 
curved waveguide arrays [2,3]. The dynamic localization 
has been confirmed experimentally in [4,5]. The same dy- 
namic control has been applied in Bose- Einstein conden- 
sates to suppress quantum tunneling of particles trapped 
in a potential well by shaking the potential [6-8] (see also 
the review [9] ) . Periodically curved waveguides have also 
been used to create defect- free surface modes [10, 11]. 

Note that the above results are in the linear regime. 
When nonlinearity is present, fundamental properties of 
dynamic localization will be altered. It is indeed the 
case with the so-called coherent destruction of tunnel- 
ing [12,13], i.e. dynamic localization in a finite domain. 
The parameter values for tunneling suppression, which 
are isolated points in the linear systems, become finite in- 
tervals for the nonlinear counterparts [14] (see also [15]). 
The method can still be proposed for dispersion man- 
agement of, e.g., a nonlinear atomic wavepacket [16] and 
gap-solitons [17]. 

Discrete spatial solitons occur when the nonlinear- 
ity effect balances the dispersion and diffraction of the 
system. In periodically curved optical waveguides they 
were reported experimentally in [18]. Different from self- 
trapping in straight waveguides, the solitons are obtained 
after transitional self-induced beam broadening. Nonlin- 
ear surface waves were observed experimentally in [19]. 
See, e.g., [20] for a review. 

In this letter, we show that the parametric drives de- 
crease the threshold power Pth of discrete surface solitons 
reported in [21-23]. More importantly, there are critical 
amplitudes where the thresholds vanish. Nevertheless, 
curved waveguides can be parasitic when the discrete lo- 
calisation is in resonance with the drive. Finally, we show 
that symmetric and antisymmetric off-side, i.e. intersite 
and twisted, modes switch stability at approximately the 
threshold driving amplitudes for thresholdless discrete 



surface solitons. 

The optical field a„ 
is modelled by [3] 



along the propagation direction z 
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where c > and S > are the waveguide coupling 
and the nonlinearity coefficient, respectively. The defo- 
cusing case S < can obtained using the transforma- 
tion a„ — >■ (— l)"a„. The drive is xo{z) = {nsa/h)xo{z), 
where a is the separation distance between the waveg- 
uides, Hs is the substrate refractive index, xo{z) is the 
physical modulation profile, and h is the inverse of the 
light wavenumber. Here, we take xq = —A cos{ujz). By 
proper scaling, one can take S ~ lu = 1. The propa- 
gation constant q is a control parameter related to the 
constant field power P = ^^^ [flnl^- Numerically its pres- 
ence removes the commensurability requirement between 
the oscillation frequency of the wavefield and that of the 
drive. Eq. (1) is solved for periodic orbits (see, e.g., [24] 
for the methods) . The stability of a periodic orbit is ex- 
amined by calculating its Floquet multipliers. 




(a) 
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Fig. 1. Field distributions |a„| at 2 = 50 mm for varying 
A with c = 0.1. Initially only the most left waveguide is 
excited with (a) ai(0) = 0.6, (b) ai(0) = 0.2. 

To illustrate the effect of curvature on the formation of 
discrete surface solitons, we plot in Fig. l(a,b) the output 



intensity |a„(50)| for the initial condition a„^i(0) = 0, 
ai(0) = 0.6 and 0.2, respectively. These correspond to 
the case of above and below Pth of the undriven case, 
i.e. P sa 0.33 (calculated using the method of [21,23]). 
In the first case, when A — Q& discrete soliton does not 
form because the repulsive force of the surface [23] is 
large enough to push the excitation to the right. Nev- 
ertheless, a relatively small driving amplitude A sa 0.67 
is enough to form a surface soliton. More interestingly, 
for the second case we observe the formation of discrete 
solitons at particular values of A even though the power 
is far below Pth of the undriven case, which will be ex- 
plained by studying waves of (1) that are periodic-in-z 
but localized in the transverse direction. 
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Fig. 2. (a) Amplitude profile at z = for discrete surface 
solitons centered at different sites near the surface with 
c = A = 0.1 and q = 0.5. (b) The power of discrete soli- 
tons in (a) for varying q. The solid black line corresponds 
to discrete 'bulk' solitons. Dotted curves show instability, 
(c) Pth as a function of A. (d) A surface phonobreather 
for the same parameter values and q = 1. 

We depict numerically exact discrete surface solitons 
of (1) centered at the mth waveguide in Fig. 2(a) for 
c = 0.1, A = 0.1, and q = 0.5. In the presence of para- 
metric drives, we show in Fig. 2(b) the power of discrete 
surface solitons for varying q. One can observe that as q 
decreases the power also decreases. Yet, there is a critical 
value of q where the power starts increasing. This implies 
that there is a minimum power for surface modes to ex- 
ist, as for the undriven case. As a comparison, we also 
present as solid black line the power variation for dis- 
crete 'bulk' solitons that have no threshold power, i.e. 
Pth = 0. We also calculate the stability of the modes 
where there is a pair of Floquet multipliers at -1-1 at Pth- 
The multipliers leave the unit circle as the continuation 



passes the minimum power point, i.e. the discrete surface 
solitons are unstable and represented by the dotted lines 
in Fig. 2(b). 

Studying Pth further, we found that it is a function of 
A. Performing similar continuations as in Fig. 2(b) for 
various driving amplitudes, Pth is plotted in Fig. 2(c). 
One can observe that Pth decreases with the increase of 
A. The decrease in Pth is responsible for the formation of 
surface solitons observed for nonvanishing A in Fig. 1(a). 
More importantly, in Fig. 2(c) there are critical ampli- 
tudes where Pth = 0. At these values of A, we obtain 
discrete surface solitons without any threshold power. 
The formation of discrete solitons in Fig. 1 (b) occurs at 
these particular driving amplitudes. 

Despite the constructive effect, we also observed that 
the parametric drive can be parasitic by inducing radi- 
ation from the excited site. In that case, we could not 
obtain genuinely localized waves and instead found dis- 
crete surface solitons with non-zero tails as shown, e.g., 
in Fig. 2(d). Such solutions are commonly referred to as 
phonobreathers, i.e. localised waves that are in resonance 
with a phonon (see, e.g., [24] and references therein). As 
the phonon band of (1) is the interval ±[q — 2c, q + 2c], 
it is straightforward to obtain the resonance condition 
k G [q—2c,q + 2c], for an arbitrary integer k. Thus, 
discrete solitons cannot exist within these regions. 
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Fig. 3. (a-b) The same as Fig. 1, but for initial conditions 
0'n^i,2 — 0, ai^2 — vO^ (a) and ai = ~a2 — a/I15 
(b). (c-d) The stability diagram of the intersite (c) and 
twisted (d) mode for fixed power P = 1. Shown is q 
against A. Dotted lines correspond to unstable solutions. 

We have also studied symmetric and antisymmetric 
two-excited site modes, which are referred to as (in- 
phase) intersite and (out-of-phase) twisted modes. For 
discrete 'bulk' solitons in infinite arrays it is known that 



in-phase modes are always unstable and out-of-pliase 
modes have a stability interval. Figs. 3(a-b) show the 
field distributions at z = 50 of the initial state where 
only the first two left waveguides are excited in-phase 
and out-of-phase, respectively. We observe that there are 
some intervals of A where intersite discrete surface soli- 
tons are obtained. On the other hand, there are also in- 
tervals where twisted surface solitons are not obtained, 
even though the results are less clear than the first ones. 
These indicate that periodically curved waveguides can 
support in-phase discrete solitons as well as destroy out- 
of-phase modes. 

We have also sought the exact periodic intersite and 
twisted modes of (1). The upper (dotted) curves in Fig. 
2(b) correspond to intersite surface solitons. We depict 
in Fig. 3(c-d) the stability of those modes as a func- 
tion of A. The power is fixed and the same as the initial 
conditions. We observe that the drives can indeed stabi- 
lize and destabilize unstable intersite and stable twisted 
states, respectively. The stability regions of the in-phase 
solitons are in agreement with the region for the for- 
mation of two-excited site surface solitons in Fig. 3(a) 
with the stability switchings occuring approximately at 
the critical drives for thresholdless discrete surface soli- 
tons (see Fig. 2(c)). As for twisted modes, it is rather 
only the stability window of Fig. 3(d) for ^ > 5 that is 
comparable with that of Fig. 3(b). While all the stability 
switchings arc due to exponential instability, i.e. Floquet 
multipliers leave at -1-1, we observe that the first stabil- 
ity switch in Fig. 3(d) as A increases from is due to a 
Hamiltonian Hopf bifurcation. 

Analysis of (1) is usually performed through its aver- 
aged equation [3] 



zd„ = cJo[an+i + a„_i] -t- (S |a„|^ 



qan, 



(2) 



where Jo{A) is a Bessel function of the first-kind. Thus, 
A y^ decreases the effective coupling between the 
waveguides yielding smaller Pth than the undriven case. 
Therefore, thesholdless discrete surface solitons are ex- 
pected to be achieved when Jo{A) — 0. Despite the sim- 
ilarity with the condition for dynamic localizations, the 
resonance here is related to self-trapped states, i.e. no 
diffusion and dispersion of fields. As for in-phase off-side 
modes, using (2) one would expect the stability switching 
to occur when Jq < 0, i.e. effectively the modes become 
twisted modes. Using a similar observation, the twisted 
mode should have been unstable when Jo{A) > 0. Nev- 
ertheless, (2) could not predict the first instability win- 
dow in Fig. 3(d) due to a quartet of multipliers. (2) is 
expected to be valid when e^"" rapidly oscillates, i.e. 
Auj ^ 1 when all the other parameters are of 0(1). 

To conclude, we have shown numerically that periodi- 
cally curved waveguides can control the formation as well 
as annihilation of discrete (surface or 'bulk') solitons. 
The parametric drives can also be parasitic to discrete 
solitons by creating resonances with the phonon. 
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